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We propose a new model in the teleparallel framework where we consider a scalar field nonmin- 
imally coupled to both the torsion T and a boundary term given by the divergence of the torsion 
vector B = |9 M (eT M ). This is inspired by the relation R = —T + B between the Ricci scalar of 
general relativity and the torsion of teleparallel gravity. This theory in suitable limits incorporates 
both the nonminimal coupling of a scalar field to torsion, and the nonminimal coupling of a scalar 
field to the Ricci scalar. We analyse the cosmology of such models, and we perform a dynamical 
systems analysis on the case when we have only a pure coupling to the boundary term. It is found 
that the system generically evolves to a late time accelerating attractor solution without requiring 
any fine tuning of the parameters. A dynamical crossing of the phantom barrier is also shown to be 
possible. 


I. INTRODUCTION 

The late time acceleration of the universe is one of the great problems of modern physics. Ever since this acceleration 
was suggested by Type-la supernovae surveys 00, this important result has been backed up by ever increasingly 
precise cosmological observations, including measurements of the cosmic microwave background (CMB) 00] , the 
Hubble constant , baryon acoustic oscillations @, and further measurements of Type-la supernovae [7]. Despite all 
these observations, there is a severe lack of understanding of this acceleration from a theoretical point of view. The 
standard cosmological model assumes the existence of a cosmological constant to explain this acceleration. However, 
a cosmological constant suffers from theoretical difficulties dues to its extremely small observed value compared to 
predictions from quantum field theoretical considerations, see 0 for a review of this subject. 

Teleparallcl gravity is an alternative formulation of gravity which is equivalent to general relativity. General 
relativity is a geometric theory based on the Levi-Civita connection, which possesses curvature but zero torsion. On 
the other hand teleparallcl gravity uses a different connection. There is a result originally discovered by Weitzenbock 
that it is always possible to define a connection on a space such that is globally flat, in other words it possesses zero 
curvature. This connection is called the Weitzenbock connection, and although it has a vanishing Ricci tensor it 
has non trivial torsion. This result is used to formulate an action based on a gravitational scalar called the torsion 
scalar T , which uses the Weitzenbock connection. The dynamics of this action are completely equivalent to general 
relativity, and this follows from the following result 


R = -T + B 


(1) 


where R is the Ricci Scalar, and B is a boundary term related to the divergence of the torsion. Since B is a total 
derivative, it gives no contribution to the field equations, and hence the action of teleparallel gravity is completely 
equivalent to the Einstein Hilbert action. 

There are numerous ways to modify both general relativity and teleparallel gravity with the aim of being able to 
understand the late time acceleration of the universe without the need for a cosmological constant. One approach is 
to modify the gravitational sector and one can consider both f(R) theories of gravity and f(T) theories. In general 
these give rise to different dynamics, since one can no longer write f{R) = f(T ) + a total derivative. f(R) gravity 
generically has fourth order field equations. In some sense /(T) theories are less harmful in that this modification 
only results in second order field equations, however the price you pay for this is that local Lorentz invariance is lost. 
In a recent paper 0 a more general theory was considered in the teleparallel framework called f{T , B) gravity, where 
B is a divergence of the torsion, which incorporates both f(R) and f(T) gravity as particular subcases. 

Another approach is to modify the contents of the universe and add a form of matter into the universe called dark 
energy which possesses a negative pressure. This can be achieved by adding to the matter sector a canonical scalar 
field fl2B25l , known as quintessence, a phantom scalar field [26}j31j , or a combination of both of these fields called 
quintom [32M40l. p75l]. A review of these models can be found in |4l| and [40 . 
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One can also consider a coupling between the scalar field and the gravitational sector. The standard approach is 
to consider a coupling between the scalar field and the Ricci scalar, of the form £ Rq L > 2 [43l - l45| . Such a nonminimal 
coupling has motivations from different contexts. It appears as a result of quantum corrections to the scalar field in 
curved spacetimes El Si and it is also required by renormalisation considerations [45| . It also appears in the context 
of superstring theories Sq. Such models have attempted to explain the early time inflationa ry epoc h, however the 
simple model with a quadratic scalar potential is now disfavoured by the current Planck data EUSil- 

In recent years an alternative formulation has been considered where the coupling occurs between the scalar field 
and torsion of the form £Tq i 2 SI SI SUll] . This gives rise to different dynamics and interesting phenomenology, for 
example phantom behaviour and dynamical crossing of the phantom barrier. A dynamical systems analysis of these 
models were considered in [53], and the observational constraints on such models were found in EU. Other types of 
nonminimal coupling to the torsion sector have been considered, for example in [6Q|] a coupling between torsion and 
derivatives of the scalar field are considered. 

In this paper we consider a different approach where we examine a nonminimal coupling between the scalar field 
to both the torsion scalar T and the boundary term B. We note that coupling a scalar field to a boundary term 
is not a new idea, for example a nonminimal coupling between a scalar field boundary terms such as the Gauss- 
Bonnet term and higher Lovelock polynomials have previously been considered 16 ll |62|. Our theory encompasses 
both nonminimally coupled teleparallcl gravity and nonminimally coupled general relativity in suitable limits. 

We then consider in detail the dynamics of the model when we have only a coupling between the boundary term 
B and the scalar field. We use dynamical systems techniques to examine the global dynamics of the system. These 
techniques are widely used in cosmological applications since they allow one to study all possible evolutional paths 
when there is no possibility of finding an exact solution. 

This paper is organised as follows. In Section [TT] we introduce our conventions and the teleparallel equivalent of 
general relativity. In Section [IIII we briefly review previous models involving nonminimally coupled scalar fields, and 
introduce the Lagrangian and field equations of our model. In Section IIVI we derive the equations determining the 
background cosmology. Finally in Section [V] we perform a dynamical systems analysis in the case where we have 
purely a coupling between the scalar field and the boundary term. 


II. TELEPARALLEL EQUIVALENT OF GENERAL RELATIVITY 


In this section we introduce the teleparallel equivalent of general relativity. The dynamical variable in this theory 
are given by the tetrad e“, where Greek indices are spacetime indices, and latin indices are tangent space indices. 
The metric g^ is given in terms of the tetrad as 

9/iv = ( 2 ) 

where r) a b is the Minkowski metric. We also introduce the inverse tetrad E%, such that 


E 11 e n = 5 n 


and E^e™ 


= 6 v 


( 3 ) 


The quantity e is defined to be the determinant of the tetrad e“, and is equivalent to the volume element of the 
metric, e = y/—g. 

In what follows we will obey the conventions outlined in 0. In general relativity the Levi-Civita connection is 
used, however in contrast, in teleparallel gravity the Weitzenbock connection a „ is chosen, and is given by the 
following derivatives of the tetrad 


W„\ = W (4) 

This connection possesses torsion but the Ricci tensor of this connection vanishes, and hence possesses zero curvature. 
The torsion tensor is simply the antisymmetric part of the Weitzenbock connection 

T a [us = W/„ - - d v e%. (5) 

The tensor X^, which we call the torsion vector, is defined as the unique nontrivial contraction of the torsion tensor 

T m = TV (6) 

Now the field equations of teleparallel gravity follow from varying the following Lagrangian density with respect to 
the tetrad 

r _ ^ oabcrri 

Lt ~ * ±abc ' 


( 7 ) 
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where the tensor S is defined as the following 

gabc — _ rpbac rjicab ^ _j_ _!_ (^ q acr ]^b T] ab T C ^ ^8) 

This combination of S abc and T a bc is usually denoted by T, which we call the torsion scalar 

T = S abc T abc . (9) 

Now in order to show that teleparallel gravity is equivalent to general relativity, we must express quantities using 
the Weitzenbock connection in terms of quantities involving the torsion free Levi-Civita connection. One can write 
the Levi-Civita connection °T in terms of the Weitzenbock connection as so 

° r A P = (10) 

where here K is called the contortion tensor and it is defined as 

2K li x v = T x llv -T vtt x +T tl x v . (11) 

This contortion tensor is antisymmetric in its last two indices. Now expressing the Ricci scalar of the Levi-Civita 
connection in terms of the Weitzenbock connection using COD, the following relation can be found 

R = -T + ^(eT"). (12) 

As the difference between the Ricci scalar and the torsion scalar is simply a total derivative, the action 0 gives 
rise to the same dynamics as the Einstein Hilbert action. This shows that teleparallel gravity is indeed equivalent to 
general relativity. Defining the boundary quantity 


B = ^(eT'Q (13) 

one then has simply the relation R = —T + B. In Q this led us to consider the modification f(T,B) gravity. In this 
work we will follow a different approach to include the term B nontrivially in the dynamics, by coupling it to a scalar 
field. We also note that one can write B in terms of a Levi-Civita covariant derivative simply as B = 2 


III. NON-MINIM ALLY COUPLED SCALAR FIELDS 

The first approach to nonminimally coupling a scalar field to the gravitational sector is to consider a coupling to 
the Ricci tensor as follows 


S = 


7^ + + m 2 ) - + L n 


\J—g d A x . 


(14) 


This approach was originally considered in the context of Brans-Dicke theories motivated by introducing a variable 
gravitational constant depending on a scalar field. In this notation the effective gravitational constant can be written 
as 


Geff 


G 

1 + K 2 t;(j) 2 ’ 


(15) 


and one can redefine such a scalar field so that it coincides with the standard Brans-Dicke field. 

Non-minimally coupled quintessence corresponds to taking £ = 0 in the above Lagrangian. For a review of 
quintessence type models, see [4l|. Quintessence alone can give rise to many interesting features from late time 
accelerated expansion of the universe to inflation [I3-G5- However simple models of scalar field inflation are becom¬ 
ing disfavoured by the latest Planck data. Another issue with a simple quintessence approach is that the effective 
equation of state must always satisfy w e ff > — 1 and require a very flat fine tuned potential in order to explain current 
cosmological observations. 

When £ ^ 0, the nonminimal coupling can be transformed to a minimal coupling via a conformal transformation 
from the Jordan frame to the Einstein frame. Such a transformation reduces the system to a quintessence model 
with a coupling between the scalar field and dark matter. For a constant exponential potential, the system can then 
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be written as a two dimensional dynamical system. Physical quantities in this frame can then be transformed back 
into physical quantities in the Jordan frame. For a review of the dynamics of these models, see chapter 9 of [4l| and 
references therein. Alternatively, one can work directly in the Jordan frame; a dynamical systems analysis for various 
potentials have been considered by various authors, see [fijl Icg 1 | and references within. 

An alternative approach has been to consider a scalar field nonminimally coupled to torsion, giving rise to teleparallcl 
dark energy theories [52|. The following action is considered 


S = 


+ l 2 ) -VW + L„ 


edrx. 


(16) 


This gives rise to different dynamics to the case of the nonminimal coupling to the Ricci scalar. Of course with a 
minimal coupling, setting £ = 0, the two theories again become equivalent due to the teleparallel equivalence. This 
theory again has a richer structure than simple standard quintessence behaviour, with both phantom and quintessence 
type dynamics possible, along with dynamical crossing of the phantom barrier. 

The equivalence between general relativity and teleparallcl gravity breaks down as soon as one nonminimally couples 
a scalar field, the field equations result in different dynamics. In this paper, we consider a more general action [791 ]. 
with the aim of unifying both of the previous considered approaches 


S = 



2^ + - X B4> 2 ) ~ 


; d 4 a 


(17) 


When one sets X = — £ one will recover an action which is equivalent to (HU), and when one sets x = 0 the action (fl6l) 
is recovered. A particularly interesting subcase will be when £ = 0, corresponding to a pure coupling between the 
boundary term. Such a coupling has not been previously studied in the literature. One could in principle choose a 
more general coupling r/((f>)B between the potential and the boundary, however the choice ry(0) = </> 2 ensures that the 
constant X is dimensionless. 

We now derive the field equations of the action m- Varying the action with respect to the tetrad field yields the 
following field equations 




e-%(eS a n - E x T p pX S p vp - -E v a T 


-E v n 


-d^d^ - v(4) 


+E»8rw^ - 4(£ + x)E p a Sr - x E v a U(d> 2 ) - E^V^cf) 


_ rpV 


(18) 


where □ = and V is the covariant derivative with respect to the Levi-Civita connection. Here T" is the 

standard energy momentum tensor derived from varying the matter sector, and is not to be confused with torsion. 

The term in the field equations proportional to S p ^ v is Lorentz violating [Tlj . This term vanishes only in the case 
when we have £ = — x, which corresponds to the case where the nonminimal coupling reduces to that of a coupling 
to only the Ricci scalar. This result is along similar lines to that obtained in Q, where the modification f(T, B) was 
considered, where it was found that the field equations are invariant under local Lorentz transformations if and only 
if f(T,B) = f(-T + B) = f(R). 

Now it can be shown that the Einstein tensor of the Levi-Civita connection can be related to the torsion sector via 
the relation 


G a v = - 


2e~ 1 d p {eS a n - 2 E^T p pX S p ^ - -E"T 


This means we can write the field equations in a covariant form as follows 


— + 2 £d> ) Gpv — g p , 


lv x d>v x d>-v(d>) 


- 


-4(£ + x)V u<j>d\d> - X 9nv u {<t> 2 ) - V M V„(<)> 2 ) 


= T„ 


(19) 


( 20 ) 


It is readily seen that this equation reduces to the correct field equation for a nonminimal coupling of a scalar field 
to the Ricci scalar when one takes X = 

Finally we have the modified scalar field equation. This is obtained by varying the action with respect to the scalar 
field, yielding the following Klein-Gordon equation 


Ud>+V\d>) = (.iT + x B)d>- 


( 21 ) 


















5 


IV. COSMOLOGY 


In this section we will derive the background equations for the cosmology of the above models. We will consider 
the standard spatially flat Friedmann-Robertson-Walker (FRW) tetrad e“ = diag(l , a(t), a(t), a(t)) corresponding to 
a spatially flat FRW metric 

ds 2 = dt 2 — a(t) 2 (dx 2 + dy 2 + dz 2 ), (22) 

where a{t ) is the scale factor. We will assume the energy momentum tensor of the matter sector is standard barotropic 
matter given by an isotropic perfect fluid 

T? = diag(p, -p, -p, -p), (23) 

where p is the matter energy density and p is the pressure. We will also assume all dynamical quantities, including 
the scalar field (j>, are homogeneous, depending only the time t. 

Inserting this FRW tetrad into the field equations (fl8l) gives us the following Friedmann equations 

3 H 2 = k 2 (p + pt ), (24) 

3H 2 + 2H = —At 2 (p + p^ . (25) 

Here H is the Hubble parameter H = ^ and we have defined the energy density and pressure of the scalar held as 
follows 

Pi = + V{4>) - 3 tH 2 <\> 2 + 6 X H<W, (26) 

Pi = hi - 4 X h 2 - V(0) + 2ff#(2£ + 3 X ) + 3 H 2 <t> 2 (t + 8 X 2 ) + 2 + 6 X 2 ) + 2 xWtt) + 12+ x)- 

(27) 

And the Klein-Gordon equation (1211) reduces to 

ij> + 3 Hj> + 6 {^H 2 + x(3 H 2 + H))<j> + V'(<f>) = 0. (28) 


In the above derivations we have used that the torsion scalar and boundary term can be written as 

T = —6H 2 , B = —18H 2 — 6H. 


(29) 


We see that the coupling to the torsion scalar introduces a term proportional to (f) 2 in the first Friedmann equa¬ 
tion (12411 , whereas the boundary term involves the addition of a <f)(f> term. One can define the equation of state of the 
dark energy or scalar held as the following ratio of the scalar held pressure and energy density 


W4, = 


P± 

Pi 


And also we dehne the total or effective equation of state as 

P + P<t> 
p + Pi 




Analogously to the standard matter energy density 


11 rri — 


AC 2 p 

3 H 2 ’ 


we will dehne the density parameter of the dark energy or scalar held as 


3 H 2 ' 


so that the relation 1 = Q. m + holds. 

In this model matter obeys the standard conservation equations 

p + 3H(p + p) = 0, 


and hence also 


(30) 

(31) 

(32) 

(33) 

(34) 


Pi + 3 H(pi + pi) — — 0 


( 35 ) 





6 


V. NONMINIMAL COUPLING PURELY TO THE BOUNDARY TERM 


In this section we analyse in detail using dynamical systems techniques the case where £ = 0, where we purely have 
a coupling of the scalar field to the boundary term. In this case the quantity Q$ is always positive, and so there is an 
energy flow from dark matter to dark energy. 

Let us introduce the dimensionless variables 


n 2 p 2 K 2 (f > 2 

- ’ X - 6 H 2 


k 2 V 


z = 2y/6nx4> ■ 


which straightforwardly generalise the normalised variables used to analyse standard quintessence 
Friedman equation (12411 written in these variables is simply the constraint 

i 2 i 2 i 2 i 

1 = a + x + y + xz, 


(36) 
The first 

(37) 


which will define the boundary of our phase space. 

The phase space will be three dimensional, and we choose to work with the variables x, y , and z. Since the energy 
density of matter is non-negative, the relation 

x 2 + y 2 + xz < 1 (38) 

must be satisfied. As in standard quintessence, due to symmetries of the system, we can assume without loss of 
generality that our potential is positive and so we only need to consider y > 0. There is no restriction on the sign of 
x or z, since </> can be positive or negative, as can y. This means that generically our phase space is non-compact, 
except in the case of a minimal coupling with y = 0. If we further introduce the variables x = u + v and z = —2v 
then the phase space becomes 


2 2 , 2^i 

u — v +y < 1, 


(39) 


which we see is simply hyperbolic space H 2 . 

In what follows we will assume that the pressure is linearly related to the energy density via the standard equation 
of state, p = wp, with w the constant matter equation of state parameter, w is physically constrained to lie between 
w = 0 corresponding to (dark) matter and w = 1/3 corresponding to radiation. 

We define the quantity N = In a and denote derivatives with respect to A by a prime 


dx 1 dx 
dN H dt 


(40) 


Now in these variables, the equations of motion can be written as the following autonomous system of first order 
differential equations 


V6y 2 X(xz — 2) + 3(2a; + z) (a: 2 (4y + w — 1) + x(w 

(z 2 + 4) 

y (V&X (x {z 2 + 4) + 2 y 2 z) + 4 (3 x 2 (w — 1) + 3 x(w 

2 (z 2 + 4) 


z' = 12yar. 


1 )z + y 2 (w + 1) — w + l) 


1 )z + 3 y 2 (w + 1)) — 6(2w; + 2 + z 2 )) 


(41) 

(42) 

(43) 


Here we have defined the quantity A by 


V'M _ 

kV (</>) 


(44) 


In order to close the system we will have to specify a form of A. For the autonomous system to remain three 
dimensional, one needs to choose a form of the potential such that A can be written in terms of the variables x, y and 
z. In this work we will assume the potential V to have an exponential form of the kind 


= Ho e~ x ^ , 


(45) 


with Vo a constant. This ensures that A is simply a constant. However this is not the only choice that will give a 
closed three dimensional system. One could also consider a power law potential of the form 


V(0) 


M “+ 4 

(f) a 


(46) 
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where a is a constant and M is a positive constant with the units of mass. In the context of late time acceleration of 
the universe, a is usually taken to be positive, however when considering inflation negative a is often considered. In 
either case, this would allow one to write A in terms of z as 


2 V / 6 «X 

A — -. 

Z 

however we will leave the analysis of such a potential for future work. 


(47) 


A. Finite critical points 

Let us denote the above autonomous system ©-(USD as 

Xi = fi(x,y,z), Xi = (x,y,z). (48) 

Critical or fixed points of (H51) correspond to (a;*,?/*,£*) that are solutions to all three equations fi(x*, y*, 2 *) = 0. We 
note that the right hand side of the dynamical system ©-(© suffers from no singularities since the denominator is 
always well defined. 

We can rewrite physical quantities in terms of the new variables x, y and 2 . We find for the effective equation of 
state that 


12x 2 (4% + w — 1) + 12(«; 

Wef? =- 


1 )xz + 2 y 2 (\/6A z + 6w + 6) 
3 (z 2 + 4) 


3 ( z 2 + 4 w) 


(49) 


Now at a fixed point, one can integrate the second Friedmann equation (1251) explicitly to find a, and it is found 

aoc (i-i 0 ) 2/[3(1+tUeff)] . (50) 

This means that the universe expansion will be accelerating if w e g < —1/3. We can also express the energy density 
of the matter and scalar field in terms of x and y 

f l m = 1 — x 2 — y 2 — xz, (51) 

= x 2 + y 2 + xz. (52) 


Point 

X 

y 

Z 

Existence 

O 

0 

0 

0 

VA, x 

A± 

±i 

0 

0 

x = o 

B 

/s"(l+uO 

V 2 A 

ri \j (i+u>)(i-™) 

V 2 A 

0 

X = 0 and A 2 > 3(1 + w) 

C 

A 

V6 


0 

X = 0 and A 2 < 6 

D 

0 

1 

[2 \ 

V 3 A 

VA, x 

E 

0 

1 

0 

A = 0 and x/0 


TABLE I. Critical points of the autonomous system ©-(m, along with the conditions for existence of the point. Points A±, 
B and C exist only when x = 0, corresponding to standard quintessence. 

A list of the critical points, along with their condition for existence is displayed in Tab. Q] We emphasise that the 
points A±, B and C appear only in the case when y = 0, and are simply the standard critical points of quintessence 
models with an exponential potential. However, they are included in the table for completeness. In the generic case 
when y / 0 and A ^ 0 the system has only two critical points given by O and D. This is because the third equation 
of our dynamical system (l43l) forces x = 0, and thus this severely restricts the possible solution set of the remaining 
two equations. In the case of a nonminimal coupling to the torsion T there are in general three finite critical points, 
in this case the number is reduced to only two. 




















In order to determine the linear stability of these points, one must analyse the following Jacobian matrix of partial 
derivatives 


J = 


dfi{x,y,z) 

dxj 


i,j = 1,2,3 


(53) 


evaluated at each of the critical points. Information about the stability of each of the points is then contained in the 
eigenvalues of this matrix. The point will be stable if all three eigenvalues have negative real part, unstable if all three 
eigenvalues have positive real part, and will be a saddle point if there are both positive and negative eigenvalues. If 
one of the eigenvalues is zero then the critical point is said to be non-hyperbolic, and then one must go beyond a 
simple linear stability analysis treatment, for a discussion of how to do this, see for instance [73]. However, all of our 
critical points are hyperbolic so such a treatment will not be needed here. 


**Point 

Weff 

Acceleration 

Eigenvalues 

Stability 

O 

0 

No 

|,-I (VI - 16x + 1) , I (VI - 16x - 1) 

Saddle node 

A.- 

1 

No 

3 ’ 3 + \/i A 

Unstable node: A > — V6 

Saddle node: otherwise 

A+ 

1 

No 

3, 3 - J | A 

Unstable node: A < VE 

Saddle node: otherwise 

B 

0 

No 

3 3 V 24—7A 2 3 3V24-7A2 

4 ' 4A ’ 4 1 4A 

Stable node: 3 < A 2 < 24/7 

Stable spiral: A 2 > 24/7 

C 

A 2 —3 

3 

A 2 < 2 

A 2 - 3, | (A 2 - 6) 

Stable node: A 2 < 3 & A/3 > (A 2 — 3) 

Saddle node: /3 < (A 2 — 3)/A 

D 

-l 

Yes 

o 3 / ! , V* 2 -48 x+6 \ 

d ’ 2 \ i+ V* 2 +6 ) ’ 

3 ( i V a2 —48x+6 N \ 

2 V V* 2 +6 ) 

Stable spiral: 48y > A 2 + 6 

Saddle point: 0 < 48y < A 2 + 6 

Stable point: x < 0 

E 

-l 

Yes 

-3, -| (Vi^8x + 1) 
§(Vr^-i) 

Stable spiral: x > 1/8 

Stable node: 0 < x < 1/8 

Saddle point: x < 0 


TABLE II. Stability and eigenvalues of the critical points of the dynamical system (|4T1)-(|43l) along with the effective equation 
of state assuming a (dark) matter equation of state w = 0. 

The eigenvalues and stability properties of the critical points are displayed in Tab. [TT] The points A±, B and C 
exist only for x = 0, in which case z = 0 identically and the system reduces to a two dimensional system. Hence only 
the eigenvalues of the reduced two dimensional system are displayed. The eigenvalues of point D mean that there are 
three three separate cases: for 0 < 48y < A 2 + 6 all three eigenvalues are negative and so the point is a stable node, 
if 48y > A 2 + 6 there is one negative and two imaginary eigenvalues (with negative real part) and so the point is a 
stable spiral, whereas for x < 0 there is always at least one negative and one positive eigenvalue and hence the point 
is a saddle. 


B. Critical points at infinity 

Now since the phase space is not compact, we must carefully check whether there are critical points at infinity. An 
analysis of the dynamics at infinity is important to understand the global stability of the system. There are a number 
of approaches one can take to doing this. One approach is to use projective coordinates, a technique recently used 
in [631 ]. One could also use a function such as arctan to compactify the variables, see for example [68j . 































9 


In this work we follow the approach used in [53[, and use Poincare variables to compactify the phase space 
begin by introducing the compactified coordinates x r , y r and z r like so 


. We 


x r 


x _ V 

%/TT^’ Vr - VTT^ 2 ’ 


Z 

\/l + r 2 ’ 


(54) 


where r 2 = x 2 + y 2 + z 2 . We also define the quantity p = ,T _^ , so that x 2 + y 2 + z 2 = p 2 . This means the dynamics 
at infinity will now be captured by taking the limit p —> 1. We then make a further coordinate transformation, 
transforming the Poincare variables into spherical polar coordinates as so 


x r = pcosOsinip, z r = p sin 9 sin ip, y r = p cos p, (55) 

where the variables lie in the range p £ [0,1], 9 £ [0,27t] and since we are restricting ourselves to y > 0 the angle tp 
lies in the restricted range <£€[—§,§]. 

The Friedmann constraint (1381) can be written in Poincare variables as 

2x 2 + 2y 2 + z 2 + x r z r < 1, (56) 


and hence the phase space will be the intersection of this ellipsoid with the Poincare sphere 

x 2 + y 2 + z 2 < 1. (57) 

Transforming our dynamical system (141l) - (143|) into these new variables, in the limit p — > 1 we find the following 

p' = 0, (58) 

-y/l — p 2 8' = \/&\ cos 9 cos p cot tp, (59) 

\Jl — p 2 p' = \J ^ A cos p (2 sin 9 cos 2 tp + cos 0 sin 2 tp} , (60) 

and hence the angular part of the equation decouples. Setting the right hand side of these equations equal to zero, 
we find that we must have cos tp = 0, and hence the critical points are those at infinity which obey 


x r = ± COS0, 
z r = ± sind, 

Vr = 0. (61) 


Now we can use this to find an equation for 6'. We go back to the equations derived for x' r ,y' r and z' r and insert 
the ansatz m and then use the chain rule to find an expression for 9' 

9' = 24\ cos0cot0(sind + cos 9) — ^(2cot0 + l)(sin(20) + cos(26() + 1). (62) 

Setting the right hand side of this equal to zero gives the critical points at infinity. 

In the Poincare variables we have that the dark energy density parameter is given by 




x 2 + y 2 + x r z r 
1 — a: 2 — y 2 — z 2 


(63) 


At the critical points of (1621) . this dark energy density parameter is divergent. Similarly, since the relation 1 = f2 m + f20 
the matter energy density will also be divergent at this point, and so these points are not of physical interest jUJ, they 
are only of mathematical interest. Thus the model has no physically relevant critical points at infinity. Nonetheless, 
since our phase space is non-compact, and the boundary of our phase space is an abnormal shape, we will continue 
using the Poincare variables when plotting the trajectories. 


C. Cosmological implications 

In this section we will discuss the dynamics of the above system. We will first briefly review the dynamics of 
minimally coupled quintessence and nonminimally coupled teleparallel energy in order to compare our results. 
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FIG. 1. Phase space showing trajectories of standard quintessence models, for the particular parameter choice w = 0, A = 1 
(with x = 0). The point C is the late time accelerating attractor, with the shaded region indicating the region of acceleration. 


The minimally coupled quintessence is achieved by setting % = 0 in the above model. The dynamical system 
is then two dimensional and the phase space is simply the upper half unit disc. This scenario has potentially five 
critical points, O, A ±, B and C. Point O is the matter dominated point, and it is always a saddle point, as required 
by cosmological observation. The points A± are both unphysical stiff matter points with effective equation of state 
Wed = +1. These points can either be saddle points or early time attractors depending on the value of the parameter 
A. The point B is of particular interest cosmologically since it is a scaling solution, the effective equation of state 
exactly mirrors the matter equation of state w e g = w , despite the fact the energy density of the scalar field does not 
vanish at this point. 

The point C is entirely scalar field dominated, and exists only for A 2 < 6 . For a suitably flat potential A 2 < 2, this 
point can describe an accelerating universe, however the effective equation of state is bounded below by — 1 , and thus 
crossing into the phantom regime is not possible, and A has to approach zero for w e g to approach —1. For A 2 < 3 the 
late time attractor is given by this point C, whereas for A 2 > 3 point B is the global attractor. Hence it is possible in 
these models to achieve a late time accelerating attracting solution. A typical plot of the two dimensional phase space 
of these quintessence models is displayed in Fig. [Q where the parameter choice A = 1 is chosen. Here we see that 
many trajectories pass close to the matter dominated origin before passing through the shaded acceleration region 
before ending at the late time attractor at point C. 

Teleparallel dark energy, on the other hand, means we must take x = 0 but restore £ into our action (TT71) . The 
phase space analysis was first explored in [Hsj, and a further analysis, where critical points at infinity were taken into 
account was performed in [53| , and it is this analysis we review here. For generic £ and A, the system has three finite 
critical points, the origin O, along with two further points F and G (following the same point labels as in (53|). Along 
with these finite points, there are a further four critical points at infinity, denoted by K± and L±. 

The critical points at infinity are either saddle points or unstable, however the points K± can give rise to an 
accelerating universe, and thus these models are able to undergo transient inflationary periods. The points F and G 
both describe dark energy dominated points, with effective equation of state u> e ff = — 1. The point G is only a saddle 
node, however the point D is stable for A 2 < £, and thus for a large class of parameters these models have a late time 
accelerating attractor, without requiring any fine tuning. One can also get w e ff = — 1 while still having an arbitrary 
value of the potential, which is a significant advantage over standard quintessence, which requires a very small A to 
achieve such an acceleration. 

Now let us explore the features of our model. The points A±, B and C only exist in the limit x —>• 0, and of course 
the points exhibit the same behaviour as the above discussion. The point O exists also for x ^ 0, and corresponds to 
a matter dominated universe with no scalar field contribution. This point remains a saddle point for all \ and A. The 
model has a further two critical points D and E. The point E exists only when the parameter A = 0, meaning that 
the potential is simply a constant. This point is entirely dominated by the potential term, and simply corresponds to 
standard de Sitter type expansion, with the potential behaving exactly as a cosmological constant. This point also 
exists in the nonminimal teleparallel case and standard quintessence. 

The point D exists only for x 7 ^ 0, and so it is unique to this model, although its coordinates are independent of 
X- At this point the energy density from the kinetic energy of the scalar field vanishes, but it has both contributions 
from both the matter sector and the potential energy of the scalar field. It has an effective equation of state u> e ff = — 1 
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FIG. 2. Phase space showing trajectories of the dynamical system in Poincare variables when y = 1, A = 2 and w = 0. 

Point D is the global attractor. 

independent of the values of x and A. For positive y this point is always a stable spiral, independent of A and hence 
it will always describe a late time accelerating attractor solution without requiring any fine tuning. 

In Fig. [2] we display some typical trajectories in the three dimensional Poincare phase space for the particular 
parameter choice A = 2 and y = 1. The boundary of the phase space is given in Poincare coordinates by the 
intersection of (15B1) and (EH). Trajectories can pass close to the matter dominate origin O, before all trajectories end 
at the late time accelerating point D. 



-1.0 -0.5 0.0 0.5 1.0 

X 


FIG. 3. Phase space showing trajectories of the dynamical system (lllll - (l43ll projected onto the x — y plane when y = 10 3 , A = 2 
and w = 0. The points A± and B are quasi-stationary. Point D is again the global attractor. 

In Fig. [3] we display a two dimensional projection onto the x — y plane for the phase space when the parameter 
values are A = 2 and y is chosen so that it is close to zero, y = 10 -3 . In this case the critical points of standard 
quintessence, points A ±, B and C behave as quasi-stationary points. Trajectories are still attracted close to these 
points. In the plot shown, trajectories start near the early time unstable points A± and are drawn towards the quasi 
scaling solution B. However, B is no longer a true critical point and so trajectories then travel to the stable global 
attractor D, which is not present in standard quintessence. 

The dynamics of the dynamical system are less interesting from a cosmological point of view when one considers 
a negative coupling y. In this case the point D is no longer a global attractor, and trajectories are instead drawn 
towards the unphysical critical points at infinity. Such a scenario is displayed in Fig [T] The points lie on the 
boundary of the Poincare sphere, and these are the mathematical critical points at infinity corresponding to the fixed 
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FIG. 4. Phase space showing trajectories of the dynamical system 63- 651 when x — —10 3 , A = 2 and w = 0. Trajectories 
end at unphysical critical points lying at infinity. 


point of equation (1621) . Trajectories move towards the quasi-critical point B before ending at one of these points at 
infinity, but as already noted above they are unphysical since both fl m and IU are divergent at these points. Such 
models are therefore not physically viable. 


VI. DISCUSSION 

In this work we proposed introducing a nonminimal coupling of a scalar field to both the torsion scalar T and the 
boundary term B. This model incorporates both nonminimal coupling to the Ricci scalar and a nonminimal coupling 
to the torsion scalar in suitable limits. We analysed in detail the dynamics of the background cosmology when we 
have simply a coupling to the boundary term. It is found that for a positive coupling, the system generically evolves 
to a late time dark energy dominated attractor, whose effective equation of state is exactly —1. This is independent 
of the potential, and thus requires absolutely no tuning of the potential to achieve this. 


-0.7 


Weff 


- 0.8 


-0.9 


-1 


Log(a) 

FIG. 5. Plot of ui e ff against N for a typical trajectory when the parameters A = 2, y = 1 are chosen. The dashed line indicates 
the phantom barrier. 

Moreover while the system is evolving close to this late time attractor, the phantom barrier can indeed be crossed, 
a scenario impossible without the presence of the coupling. We display a plot typical of such behaviour in Fig. [5] 
It is seen that the effective equation of state can cross the phantom barrier, and indeed cross from both directions, 
oscillating around the barrier before settling at its final late time de Sitter type expansion. 

The global dynamics of these models are simpler than the case of the nonminimally coupled torsion scenario. There 
are less critical points, and there are no longer any physical critical points at infinity. Teleparallel dark energy also 
possesses saddle points describing an accelerating universe and hence can exhibit transient periods of inflation. Such 
a scenario is not possible in our model as we have no accelerating saddle points. 
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An investigation of the full phase space when a coupling to both T and B is present could in principle be analysed 
in future. In our case however we focused only on a coupling to the boundary term here to analyse the effects of this 
new contribution. One should also investigate the role solar system and other observational constraints place on these 
nonminimal coupled models, analogous to [-54| . In this work we focused on exploring the bad ound cosmology and 
so furthermore the cosmological perturbations should also be investigated, along the lines of 
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